Long-range-interactions induced ordered structures in deposition processes 
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We present a new model of sequential adsorption in which the adsorbing particles experience 
dipolar interactions. We show that in the presence of these long-range interactions, highly ordered 
structures in the adsorbed layer may be induced at low temperatures. The new phenomenology 
manifests through significant variations of the pair correlation function and the jamming limit, with 
respect to the case of noninteracting particles. Our study could be relevant in understanding the 
adsorption of magnetic colloidal particles in presence of a magnetic field. 

PACS number(s): 68.45.Da, 81.15.-z, 82.70.Dd, 68.10.Jy 



The study of the irreversible adsorption of colloidal 
particles onto a surface has been since long time subject 
of a great deal of interest, due to its potential applications 
to physical, physico-chemical, and biological problems 
jjj. Our global understanding of the process has been 
possible through the formulation of different models, an- 
alyzed either numerically or analytically. These models 
share a sequential and irreversible nature, and differ in 
the rules by which the particles accommodate when ar- 
riving at the surface. The various rules are responsible 
for the different values of the relevant quantities describ- 
ing the adsorbed phase, such as the jamming limit, the 
pair correlation function, and the local variance of the 
number of deposited particles. In the random sequen- 
tial adsorption (RSA) model J^j^l particles are placed at 
random positions on the substrate. If an incoming parti- 
cle overlaps with a previously adsorbed one, it is rejected 
and a new one probed; otherwise it becomes irreversibly 
adsorbed. In the ballistic model (BM) [|7|-[To[ the parti- 
cles descend to the surface following straight vertical tra- 
jectories. An incoming particle that does not reach the 
substrate is allowed to roll over the previously adsorbed, 
following the steepest descent path, until it reaches a sta- 
ble position. Only particles that fail to gain the surface 
are finally rejected. 

All these models, and their subsequent extensions, 
have been mainly implemented by considering short 
range — hard core — interactions among particles. With 
the exception of the analysis of the role played by elec- 
trostatic interactions jOJ, the case of long-range inter- 
actions remains essentially unexplored. Our purpose in 
this Letter is to analyze comprehensively the influence 
that these interactions have in the kinetics of deposition 
in a simple numerical model. We will show that when 
they are taken into account, a new aspect of the problem 
emerges. The structure of the adsorbed layer changes 
considerably, giving rise in some cases to the appearance 
of a higher degree of order in the substrate. 



To illustrate this point, we will focus on the case of 
anisotropic dipolar interactions, and present numerical 
simulations of the adsorption process on a line (1 + Id) 
and on a plane (1 + 2d). In our simulations we consider 
the adsorption of spherical magnetic particles of diame- 
ter a and magnetic moment fl — fiu, with \i being the 
magnetic moment strength and u a unit vector oriented 
along its direction. Orientation of dipoles is restricted to 
the space in which they diffuse: IR 2 in (1 + Id) and K 3 in 
(1 + 2c?) . The dipolar interaction between two particles 
i and j, located at positions f\ and fj, respectively, is 
Uij = fi 2 Eij/a 3 , with Eij being the dimcnsionless energy 

E i:j = a 3 {u, ■ Uj - 3(u t ■ fij ) (Uj ■ nj ) j r\ } / r 3 - , ( 1 ) 

where fij — fi — fj. Particles are released at random posi- 
tions over a certain initial height Zj n , having assigned an 
initial orientation Each particle undergoes a random 
walk until it either becomes adsorbed or moves away from 
the substrate a distance greater than z ou t, in this case, 
the particle is removed and a new one released. The ef- 
fect of the interactions is taken into account by means of 
a Metropolis algorithm jT^jlJ] . Suppose there are N —1 
adsorbed particles, located at points r*j, i = 1, . . . , N — 1. 
At some time t an incoming dipole occupies the position 
f)v and has a total energy E — YliLi ^Ni- At time 
t + 1 we compute a new position r^; the particle arrives 
there by means of jump of length 5 in a direction chosen 
at random. The movement to rjtj is performed rigidly, 
without changing Ujy. The energy experienced in the 
new position is E* , and the total change in the energy 
due to the movement is AE = E* - E. If AE < 0, 
the movement is accepted; otherwise, it is accepted with 
probability p = exp(—AE/T r ), where T r = a 3 k^T/ is 
a reduced temperature. Note that with this procedure 
we do not take into account gravitational forces. In or- 
der to speed up the algorithm, when an incoming particle 
reaches a position very close to the substrate (less than 
one diameter) it is attached according to the BM rules. 
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FIG. 1. Pair correlation function for dipolar adsorption, in 
the free limit T r — » oo, and for models REL and PAR at 
T r = 0. 
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FIG. 2. Chain-length density function for dipolar adsorp- 
tion, in the free limit T r — > oo, and for models REL and PAR 
at T r = 0. 



This procedure indeed assumes the presence of a strong 
short-ranged attraction between the substrate and the 
particles, which overcomes dipolar interactions at short 
distances. 

The description of the algorithm is completed by pre- 
scribing the dynamics of the dipoles. We have consid- 
ered two possibilities under scrutiny, defining two differ- 
ent models: (i) Relaxed dipoles (model REL). The ini- 
tial orientation vff of the dipoles is assigned at random. 
After every accepted movement, the moment of the ran- 
dom walker is oriented along the direction of the total 
magnetic field on its position. After adsorption, a newly 
attached dipole undergoes a last relaxation, and its di- 
rection does not change any more. There are two natural 
choices for the orientation of the first adsorbed dipole: ei- 
ther select it randomly or parallel to the substrate. Both 
cases provide comparable results, within the error bars. 
The results reported here correspond to a random initial 
dipole. (ii) Parallel dipoles (model PAR). All dipoles are 
parallel to the substrate, pointing rigidly towards a pre- 



determined fixed direction. Both models are actually the 
limits of a process of adsorption of relaxing dipoles in the 
presence of a magnetic field. Model REL is recovered in 
the limit of zero field, whereas model PAR corresponds 
to a very strong field. 

In (1 + Id) we have simulated the adsorption of dipoles 
at different temperatures onto a substrate of length L. 
The results reported here are mainly for L — 150a; larger 
system sizes provided systematically equivalent results. 
We have selected the parameters 5 — a/2, Zi n = 15a, 
and z ou t — 20a. These are conservative estimates; since 
dipolar interactions decay in average as r -3 , the energy 
diminishes more than 3 orders of magnitude from the 
neighborhood of the substrate to the launching height 
z in . Higher values of z in and z out were tested, provid- 
ing comparable results but dramatically increasing the 
computation time. 

The presence of dipolar interactions is expected to in- 
duce changes in the structure of the adsorbed phase. We 
have investigated this point in (1 + Id) by computing 
the density at jamming 9^, the pair correlation function 
g(r), and a new quantity, the chain-length density func- 
tion P(n). The BM rule adopted for the final allocation 
of the dipoles allows an incoming particle to roll over 
another, adsorbing therefore close together and forming 
chains — sets of contiguous particles — in (1 + Id) 15 1. We 
define P(n) as the density of chains of length n, per unit 
length of substrate. This function is related to the jam- 
ming limit through 9^ = J2^Li nPin). 

In order to first check our algorithm, we have gener- 
ated 1000 jammed configurations on a linear substrate of 
L = 750a, in the limit T r — > oo (free BM). We recover a 
jamming limit 9 B D M — 0.8079 ±0.0004, in excellent agree- 
ment with theoretical predictions 9^ w 0.80865 ||. The 
REL and PAR models render in the most extreme case 
of T r — (limit of zero temperature or infinitely strong 
interactions) the jamming limits 6^ EL = 0.855 ± 0.002 
and 9^ AR = 0.752 ± 0.002, respectively. The difference 
with the BM value is rather notable, larger than a 5.7% 
and 7.0%, respectively. This feature should be contrasted 
with the small variations in the jamming limit due to the 
consideration of hydrodynamic interactions in the pro- 
cess (around 1.4%) @. 

Figure ^ depicts the pair correlation function g(r) for 
BM, and both models REL and PAR at T r = 0. Model 
REL enhances greatly the maxima of g(r), which occur 
at distances r p = ap, for integer p (that is, correspond- 
ing to interparticle distances equal to a multiple of the 
diameter). We find in this case that the peaks decay 
exponentially, defining a correlation length of the form 
g REL (r p ) ~ exp(-7p/f). We estimate £ = 4.73±0.04. In 
model PAR, on the other hand, correlations are strongly 
suppressed; the peaks in the function are almost invisible 
for r > 2a. 

Given the expression (Q) for the dipolar energy, it is en- 
ergetically favorable for the particles to form linear chains 
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FIG. 3. Jamming limit 9oa as a function of the reduced 
temperature T r for the models REL (circles) and PAR 
(squares). 
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FIG. 4. Average chain length (n) as a function of the re- 
duced temperature T r for the models REL (circles) and PAR 
(squares). 



of aligned dipoles. This effect is studied in Fig. |[ where 
we have plotted P(n) for BM, and both REL and PAR at 
T r = 0. We observe again an exponential decay for REL, 
P REL (n) ~ exp(— n/n), with an estimated correlation 
size n = 4.67 ± 0.05, in good agreement with the char- 
acteristic length for the correlation function, £ = 4.73. 
This exponential behaviour, on the other hand, is lost in 
BM and PAR; in those cases, the function P(n) decays 
much faster. The most characteristic property of Fig. || 
is the way in which the weight of the distribution moves 
towards larger chains in model REL (up ton — 24, which 
is almost 1/6 of the total substrate length). Large chains 
are, on the contrary, strongly suppressed in model PAR 
(only n — 5, 1/30 of the substrate length). The witdh of 
the chain-length distribution can be estimated by means 
of the average chain length (n) , 



(n) 



(2) 



We observe the values {n) BM = 1.857 ±0.003, (n) PAR = 
1.298 ± 0.007, and (n) REL = 2.52 ± 0.02. 

The next point we have analyzed is the dependence 
of the relevant quantities of the models on the temper- 
ature. In Figure |^ we have plotted the jamming limit 
0oo for models REL and PAR, as a function of the re- 
duced temperature T r . A similar plot is made in Fig. |] 
for the average chain length (n). From these two fig- 
ures we conclude that, when increasing T r , 9^ and (n) 
jump rather discontinuously towards the limit value at 
high temperature (BM). This behaviour contrasts with 
the results reported for the diffusion-limited aggregation 
(DLA) of particles with dipolar interactions |l^,|l^]. In 
those works, the fractal dimension of the agregates was 



computed as a function of a reduced temperature T r . 
It was found that the fractal dimension of dipolar DLA 
also changes between a high and a low temperature val- 
ues, with a transition that seems to be continuous, taking 
place during almost 4 orders of magnitude in T r . In the 
present case, however, we observe a transition ocurring 
in a range of about 1.5 orders of magnitud in T r . It seems 
therefore more likely the presence of a crossover in dipolar 
adsorption. The discontinuous nature of this crossover is 
smoothed by finite-size effects, but would be more evi- 
dent in larger systems. We note that the crossover would 
happen at different characteristic temperatures T c for the 
two models. In view of Figs. an d [|, we estimate the 
values T L REL « 10" 125 and T™ R n 10-° 5 . 

Finally, in order to study the effects of dimensional- 
ity, we have considered the extension of our model to 
(1 + 2d), by adsorbing dipolar particles onto a plane of 
size L x L. In the case of free BM, we recover a jamming 
limit 0oo = 0.609 ± 0.003, in good agreement with previ- 
ous simulations ]9],[To| . When introducing dipolar interac- 
tions at T r < oo, the algorithm turns out to be extremely 
time consuming, specially for model PAR. We report here 
preliminary results only for model REL at T r = 0, on sur- 
faces of a fairly small size (L = 10a). Figure || depicts 
typical configurations in free BM and REL. The arrows 
drawn correspond to the projection onto the adsorption 
plane of the originally tridimensional dipoles. We can 
observe the local effect of the interactions, ordering the 
dipoles of neighboring particles in order to minimize the 
dipolar energy (|l|). The estimated value of the jamming 
limit is in the case REL ^ = 0.622 ± 0.003. The differ- 
ence with respect to the free case is 1.8%, smaller than 
that observed in (1 + lei). This is a result to be expected, 
since increasing the dimensionality increases the possibil- 
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FIG. 5. Jammed configurations in a square of size L — 10a. 
a) Free BM. b) REL model. In this last case one can appreci- 
ate the relative tendency of neighbor dipoles to orient parallel 
to each other. 

ities for disorder in the system. Therefore, the relative 
packing of the final state is expected to be smaller than 
in (1 + Id). A similar effect was also observed in 3d sim- 
ulations of dipolar DLA |L3| . 

In summary, we have shown that the presence of long 
range dipolar interactions during the adsorption process 
of particles on a substrate induces important changes 
in the structure of the adsorbed phase, reflected in the 
jamming limit, correlation function, and average chain 
length. In (1 + Id), one of the two models considered, 
namely relaxed dipoles (REL), shows long-range order, 
manifested in the tendency of the particles to form long 
chains, increasing therefore the packing of the substrate. 
On its turn, the model of parallel dipoles (PAR) produces 
smaller chains and results in a more disordered jammed 
phase. These results can be regarded as an effect of the 
different dynamics of the dipoles. Given the form of the 
dipolar energy (|l|), the adsorbed rigid dipoles in model 
PAR tend to repel an incoming particle approaching nor- 
mally the surface; this fact induces an average repulsion 



between particles and substrate. On the other hand, the 
relaxing dynamics in model REL results in an induced at- 
traction between dipoles and surface, attraction respon- 
sible for the larger surface coverage. When decreasing 
temperature, we observe in model REL a transition from 
a disordered phase at high T r to an ordered phase emerg- 
ing at low temperatures; in model PAR the transition is 
reversed. In the (1 + 2d) geometry the effects are less 
pronounced, reflected in a slight variation of the jamming 
limit and an induced orientational order in the dipoles. 
Our findings may open further perspectives for the de- 
velopment of new models accounting for the phenomena 
we have reported and the different properties of the sub- 
strate, as well as a technical hint to achieve an increased 
efficiency of the adsorption processes, allowing the pos- 
sibility of covering with particles a larger fraction of the 
surface. 
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